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摘要 :基于 遥感 影像 的 变化 检测 是 当前 的 研究 热点 ,可 为 区 域 生态 环境 保育 、 资 源 管理 与 燃 展 规划 等 提供 决策 支撑 。 目 前 遥感 
影像 的 变化 检测 多 基于 两 个 时 相 , 不 能 充分 地 反映 植被 在 时 间 维 的 连续 变化 特征 。 通 过 引入 信息 论 多 如 出 了 利用 时 间 信 息 炳 来 
综合 表征 植被 长 时 间 序 列 的 变化 特征 。 人 研究 以 延 河 流域 为 试验 区 ,基于 MODISZNDVT 数 据 ,应 用 时 间 信 息 业 方 法 来 计算 了 
2000 一 2010 年 该 区 域 的 植被 覆盖 变化 信息 ,厘清 了 时 空 变化 特征 。 研 究 结果 表明 和 ii 10 年 延 河 流域 的 植被 覆盖 的 变化 以 增加 
为 主 , 占 流域 面积 的 80.7% ;植被 覆盖 明显 增加 的 区 域 占 流域 面积 13.99%s 主 要 分 布 在 流域 的 东北 部 和 东南 部 ;植被 覆盖 减少 的 
区 域 占 比 2.4% ,主要 分 布 在 流域 的 西部 和 西北 部 ;严重 减少 的 区 域 占 比 1.1 久 ,主要 分 布 在 流域 的 中 部 和 西南 部 ,是 需要 重点 的 
生态 恢复 与 治理 区 域 。 时 间 信 息 炉 方法 与 回归 分 析 法 相 比 ,能 够 更 为 客观 地 表征 长 时 间 序 列 植被 覆盖 的 连续 变化 强度 和 变化 
趋势 ,可 为 区 域 生 态 环 境 的 保护 和 管理 提供 更 为 科学 的 理论 依据 8 
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Abstract: Change detection based on remote sensing images is a hot research topic, which can provide decision support for 
regional ecological conservation resource management and development planning, etc. However, the current studies mostly 
focus on jio-temporal change detection, which cannot adequately reflect the continuous change characteristics in time 
dimenéion. Jathis paper, through the introduction of information theory, we propose a novel method of temporal information 
efifiopy toheharacterize the long time change features of vegetation comprehensively. Taking Yanhe watershed as the study 
aiea and” based on Normalized Difference Vegetation Index ( NDVI) data, we used this method to calculate the change 
information of vegetation coverage in Yanhe watershed from 2000 to 2010 and identified its spatio-temporal change 
characteristics. The results showed that the change level of vegetation coverage in Yanhe watershed was mostly increased in 


this period, accounted for 80.7% of the total area. Obviously increased area accounted for 13.9% and mainly distributed in 
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northeast and southeast of the watershed. Vegetation coverage decreased area accounted for 2.4%, mainly distributed in the 
western and northwestern ;Severely decreased area accounted for 1.1%, mainly distributed in the central and southwest of 
the study area, where needs strict environmental management. Comparing with the regression analysis method, temporal 
information entropy can be more objectively in reflecting the continuous change intensity and trend of vegetation coverage in 
long time series which can provide scientific theoretical references for the protection and management of regional ecological 


environment. 
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基于 遥感 的 地 物 目 标 检测 是 评价 流域 生态 环境 可 持续 发 展 的 重要 方法 ,其 检测 方法 多 梯 ) 检测 精度 仍 有 
提高 空间 。 土 地 利用 变化 ,城市 发 展 与 生态 环境 变化 等 关系 到 人 民生 活 和 人 类 社会 可 持续 发 展 的 一 系列 问 
题 ,集中 起 来 可 归结 为 地 物 目标 的 变化 检测 问题 。 遥 感 技术 因 具 有 时 效 性 强 、 履 盖 范 围 雇 ` 可 重复 观测 等 
特点 ,可 以 定量 监测 地 表 发 生变 化 的 位 置 . 过 程 和 幅度 ,而 成 为 变化 检测 最 主要 的 技术 手段 。 变 化 检测 也 
是 目前 遥感 应 用 研究 中 的 热点 “5 。 般 感 影像 变化 检测 是 从 不 同时 期 的 遥感 数据 叱 5 定量 地 分 析 和 确定 地 表 
变化 的 特征 与 过 程 !'" 。 近 些 年 来 ,相关 学 者 相 继 提 出 了 许多 基于 遥感 影像 的 变化 检测 访 法 ,并 从 不 同 的 角度 
总 结 和 归 类 。Singh 等 根据 是 否 对 遥感 影像 分 类 将 变化 检测 方法 分 为 直接 氏 较 法 和 分 类 后 比较 法 ;Rensink 
等 里 根据 研究 对 象 层次 的 不 同 将 变化 检测 方法 分 为 基于 像素 级 此 于 特征 级 、 基 于 决策 级 3 个 层次 ;Lu 等 '” 
按照 采用 数学 方法 的 不 同 将 变化 检测 技术 分 为 算术 运算 法 , 变 兹 法 分 类 比较 法 、 高 级 模型 法 、GIS 集成 法 、 视 
觉 分 析 法 和 其 他 方法 等 七 类 。 然 而 ,从 检测 变化 的 本 质 策略 下 可 将 变化 检测 方法 分 为 双 时 相 变 化 检测 和 时 间 
序列 变化 检测 两 大 类 。 前 者 是 对 两 个 时 相 遥 感 影 像 的 变化 检测 ,常用 的 方法 有 影像 代数 法 、 主 成 分 分 析 
法 、 变 化 矢量 分 析 法 纹理 特征 分 析 法 ,马尔 科 夫 随机 场 模型 法 数据 融合 法 等 。 如 ; 邓 小 炼 等 "利用 变化 矢 
量 分 析 法 检测 了 华北 平原 的 土地 利用 和 覆盖 变化 ,取得 了 较 好 的 效果 ;Tsarouchi 等 ”人 利用 主 成 分 分 析 法 对 
TM 影像 进行 了 谈 化 检测 处 理 和 检测 结果 分 析 ;Brozzone 等 将 马尔 可 夫 随 机 场 模 型 法 引入 到 变化 检测 中 ， 
提取 了 两 时 相 影 像 中 的 空间 邻 域 信息 & 结 果 显示 对 噪声 起 到 了 很 好 的 抑制 作用 ,优化 了 检测 效果 ;Camps-Valls 
等 中 将 数据 融合 法 应 用 到 多 时 相通 感 影像 的 变化 检测 中 ,并 客观 评价 了 检测 结 

长 时 序 植被 覆盖 的 连续 动态 变化 检测 是 生态 学 研究 的 重要 领域 和 全 球 变化 研究 的 热点 。 遥 感 技术 与 信 
息 论 的 发 展 为 时 间 序 列 的 变化 检测 奠定 了 理论 与 数据 基础 。 随 着 遥感 数据 的 积累 ,长 时 间 序 列 的 遥感 影像 完 
整地 记录 了 地 表 的 变化 过 程 。 然 而 ,目前 的 研究 多 为 双 时 相 变 化 检测 。 通 过 时 间 序 列 变化 检测 可 以 充分 
地 挖掘 地 物 在 时 间 维 的 变化 信息 ,能 够 更 好 地 反映 其 时 空 变化 规律 “” 。 时 间 序 列 数据 在 地 表 覆 盖 变 化 的 
研究 中 有 广泛 的 应 用 ,其 中 最 为 直接 的 是 监测 地 表 植 被 的 变化 “2 。 常 用 的 植被 覆盖 变化 检测 方法 主要 有 
回归 分 析 治 \ 主 成 分 分 析 法 、 传 里 叶 变 换 、 小 波 变换 等 。 回 归 分 析 法 是 考虑 变量 之 间 统 计 联 系 的 一 种 重要 方 
法 ,目前 在 植被 覆盖 的 长 时 序 变化 趋势 研究 中 应 用 最 为 广泛 :” ”1 。 该 方法 的 优点 是 计算 过 程 简单 清晰, 结 
果 直 观 。 但 是 要 求 研究 对 象 需 满足 正 态 分 布 , 且 对 遥感 影像 的 噪声 .误差 规避 能 力 较 差 '” ;此 外 ,该 方法 也 很 
难 客观 地 反映 出 长 时 序 植被 覆盖 的 变化 强度 特征 。 主 成 分 分 析 法 \ 傅 里 叶 变 换 及 小 波 变 换 类 似 ,是 通过 对 时 
序 遥 感 数据 进行 影像 变换 或 数学 处 理 来 增强 变化 信息 ,然后 从 变换 分 量 中 检测 变化 信息 ,此 类 方法 可 以 有 
效 地 突出 变化 信息 ,能 够 较 好 地 表达 植被 覆盖 的 宏观 变化 特征 ,但 通过 变换 处 理 得 到 的 影像 意义 不 够 明确 , 且 
对 长 时 序 植被 覆盖 在 时 间 维 的 具体 变化 特征 难以 体现 出 来 。 鉴 于 常用 分 析 方 法 在 长 时 间 序 列 植 被 覆盖 变 
化 检测 中 存在 的 不 足 ,本 文通 过 引入 信息 论 , 提 出 了 时 间 信 息 入 的 概念 及 计算 方法 ,希望 为 该 领域 的 研究 提供 
新 思维 .新 方法 。 

时 间 信 息 粹 是 检测 植被 履 盖 的 时 序 变化 特征 的 有 效 方法 。“ 烂 ”最早 起源 于 物理 学 ,用 来 反映 热力 学 系 
统 的 无 序 程度 。1948 年 , 香 浓 将 这 一 概念 引入 通信 和 领域 并 称 之 为 信息 炉 , 它 反 映 了 信息 源 的 平均 不 确定 
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性 "3 。 如 今 ,信息 粒 已 被 广泛 应 用 到 生态 学 .地 理学 .社会 科学 等 领域 “” 。 信 息 烂 在 遥感 影像 变化 检测 领 
域 中 也 有 较 多 的 应 用 ,但 多 是 对 两 个 时 相 变 化 信息 的 检测 。 如 : 李 亮 等 提出 了 利用 般 感 影像 的 像 斑 差 炉 来 
检测 两 期 地 表 覆 羡 变 化 情况 ; 任 晖 等 “提出 了 土地 利用 异动 粹 来 提取 两 期 影像 的 土地 利用 变化 信息 ,取得 了 
较 好 的 效果 , 且 该 方法 对 所 用 遥感 数据 源 的 配 准 精度 分辨 率 等 要 求 较 低 。 本 文 在 已 有 研究 成 果 的 基础 上 , 基 
于 信息 论 ,提出 了 利用 时 间 信 息 焙 方 法 来 深入 挖掘 时 间 序 列 遥 感 影像 的 植被 覆盖 变化 信息 ,综合 、 定 量化 地 表 
征 其 在 时 间 维 的 变化 规律 。 并 以 延 河 流域 为 试验 区 ,应 用 时 间 信 息 业 方 法 来 客观 反映 延 河流 域 长 时 间 序 列 杆 
被 覆盖 的 时 空 变化 特征 ,以 期 为 该 区 域 生态 环境 的 保育 和 治理 提供 决策 支撑 。 


1 试验 区 概况 与 数据 来 源 


1.1 试验 区 概况 

延 河 流域 生态 环境 问题 较为 突出 ,是 陕西 省 生态 环 
境 建 设 的 重点 区 域 。 延 河 是 黄河 的 一 级 支流 ,发源 于 榆 
林 市 靖 边 县 天 赐 湾 乡 ,由 西北 向 东南 依次 流 经 靖 边 县 、 
志 丹 县 . 安 塞 区 .宝塔 区 和 延长 县 ,最 后 在 延长 县 南河 沟 
乡 凉水 岸 附近 注入 黄河 ,干流 总 长 286.9km。 延 河流 域 
位 于 陕 北 黄土 高 原 腹 地 ,地 理 位 置 36"23' 一 37"17'N ， 
108"45' 一 110?28 中 ,流域 总 面积 约 7687km (图 1) 。 地 
势 起 伏 较 大 ,海拔 范围 475 一 1787m ,平均 海拔 1213m。 2 
截止 2010 年 底 ,流域 总 人 口 约 95 万 ,其 中 农业 人 时 占 109°00’ 109°30’ 110°00’ 110°30 福 


37°00N 


36°30" 


61% ,GDP 总 量 为 382.7 亿 元 '”。 宝 塔 区 是 我 国 著名 图 1 延 河 流域 位 置 示意 图 
的 红色 旅游 胜地 。 安 塞 区 矿产 资源 丰富 ,石油 产业 是 其 Fig.l Location of Yanhe watershed 


支柱 产业 。 该 流域 属于 典型 的 上 暧 温带 大 陆 性 半 干 旱季 
风气 候 , 冬 季 寒 冷 干燥 ,夏季 高 温 多 雨 。 从 西北 向 东南 ,降水 .温度 具有 明显 的 梯度 变化 特征 。 年 降水 量 少 且 
季节 分 配 不 均匀 ,年 均 降 雨量 为 520mm ,年 均 气温 约 8.8 一 10.2%C 。 该 区 域 的 植被 覆盖 状况 从 南 向 北 呈 现 着 较 
为 明显 的 变化 。 科 学 .客观 地 反映 其 换 时 间 范 围 植 被 覆盖 的 变化 情况 ,可 为 该 区 域 生 态 环境 的 建设 和 治理 
提供 决策 信息 。 
1.2 ”数据 来 源 

归 一 化 植被 指数 (NDVI) 是 本 研究 的 重要 数据 源 ,其 能 够 很 好 地 反映 地 表 的 植被 覆盖 情况 。NDVI 是 还 
感 卫 星 监测 植被 应 用 最 为 广泛 的 参数 中 。 研 究 中 所 采用 的 NDVI 数据 来 源 于 美国 NASA ( https://ladsweb. 
nascom.nasa.gov/data) 的 MODIS 植被 指数 产品 MOD13Q1 ,空间 分 辨 率 为 250m ,时 间 分 辨 率 为 16d。 该 数据 产 
品 统一 采用 SIN 听 影 并 已 进行 去 云 .辐射 校正 等 预 处 理 。 人 研究 所 选用 的 时 间 范 围 为 2000 年 2 月 到 2010 年 12 
月 。 利 由 MOBLAND 提供 的 MRT( MODIS Reprojection Tools ) 软件 对 该 数据 产品 进行 格式 转换 和 投影 转换 ,将 
HD 格式 转换 为 Tif 格式 ,并 将 SIN 投影 转换 为 WGCS84/Albers Equal Area Conic 投影 , 重 采样 方法 为 最 邻近 
法 s 重 采样 分 辨 率 为 250m。 然 后 利用 MVC( Maximum Value Composites ) 法 对 各 年 23 期 数据 进行 最 大 值 合成 ， 
来 进 主 步 消除 云 . 大 气 和 太阳 高 度 角 等 因素 对 遥感 影像 产生 的 影响 ' 引 ,得 到 试验 区 2000 一 2010 年 逐年 的 
NDW 年 最 大 值 数据 。 


2 ”研究 方法 


2.1 信息 炳 
信息 箭 (也 称 狭 义 信 息 和 ) 的 计算 方法 基于 离散 型 随机 变量 ,其 计算 公式 见 式 (1) 。 


HO) =- 之 po) logsp(%) a) 
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式 中 , x; 表示 随机 变量 X 的 状态 或 取 值 ; p(x,) 为 相应 的 概率 。 
连续 型 随机 变量 的 信息 箭 不 是 对 上 式 简单 地 取 极 限 值 ,而 是 以 微分 灶 的 形式 定义 的 , 见 式 (2) 。 


HOX) =- [f(s) logsf(x) dx (2) 


式 中 , f(x) 为 概率 分 布 函 数 。 

由 上 面 的 计算 公式 可 知 , 求 信息 粹 最 关键 的 是 确定 概率 密度 或 者 具体 的 概率 分 布 函数 。Vasicek 等 '“ 提 
出 了 直接 基于 观测 数据 ,而 无 需求 概率 密度 函数 或 概率 分 布 函数 来 计算 信息 炳 的 估计 值 ,并 证 明了 该 估计 值 
依 概率 收敛 到 理论 值 。 将 式 (2) 表 示 为 式 (3) 的 形式 : 


#0 = fA) tog) dr = los (G77) ) (3) 
式 中 ,Fp-1(p) 的 微分 可 由 样本 点 处 的 斜率 值 来 近似 表达 , 见 式 (4) ; 
Virm Yi-m 
h,(m,n) = ee (4) 


式 中 ,yy 夺 … 二 yy ,是 样本 点 x ,… ,x, 按照 从 小 到 大 的 顺序 排列 得 到 的 。 知 为 正 整 数 , 且 0 < m= 本 。 


进而 ,信息 箭 可 以 表达 为 式 (5) : 


一 人 
H = ] 5 
(m,n) = > og, "0 BA C5) 
Ebrahimi 等 524 对 上 式 进行 了 改进 ,并 提出 IE 见 式 66) : 
本 JS 区 a 
| H S l 6 
em si = > pe2 x mn ci Xm/n 8 
i—-1 
1 + ,] im 
L <m 
其 中 ,c=I2,m+1<i<n-m ; 了 
Yirm = Inn-m 
7 一 1 
] 十 NR-m+l1 im 


二 式 中 ;y, <y, <…<y ,是 样本 点 x ,… ,x 按照 从 小 到 大 的 顺序 排列 得 到 的 。m 为 正 整 数 , 旦 0 <m 了。 


2 
同时 ,证 明了 该 方法 依 概 率 收 敛 到 理论 值 8(X) ,并 且 与 式 (5) 相 比 有 较 小 的 偏差 和 均 方 根 误差 ”7 。 
2.2 ”时 间 信 息 信 
在 借鉴 信息 炉 概 念 的 基础 上 ,并 结合 遥感 数据 源 的 特点 ,本 文 提 出 了 新 的 概念 和 方法 ,即时 间 信息 炉 .时 
间 序 列 信 息 炉 充分 挖掘 植被 覆盖 在 时 间 维 的 变化 特征 ,表征 其 在 某 一 时 期 的 变化 强度 和 变化 趋势 信息 。 
时 间 信 息 精 (iemporal information entropy) 以 公式 (6) 为 基础 ,用 来 反映 植被 覆盖 在 时 间 维 的 变化 强度 特 


征 。 ) 冰 用 诗 应 天 芭 不 同 数据 源 的 特点 ， 0 的 NDVI 为 例 ,其 具体 计算 如 式 (7) : 
im Yi-m 
be xm 2 (7) 
i—-1 . 
1 十 ,1 i<m 
c=12 ,m+1l<i<n-m ; 后 Nm 
) yy = in-m 
1 wm 
式 中 ， yy, Ey, ,是 样本 点 NN, 按照 从 小 到 大 的 顺序 排列 得 到 后 得 到 的 ( 即 表示 像 元 第 ;年 的 
NDVI 值 )。 
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下 面 对 计 算 公 式 中 的 参数 进行 说 明 : 

m 是 “时 间 频 率 ” 因子 ,通过 选取 不 同 的 m 值 (m 为 不 超过 n/2 的 正 整数 ) 可 以 反映 人 研究 对 象 在 不 同时 间 
尺度 上 的 变化 特征 。 以 本 研究 为 例 ,由 于 试验 区 的 植被 覆盖 情况 是 用 每 年 NDVI 数据 的 最 大 值 来 反映 的 , 故 
研究 对 象 变化 的 时 间 尺 度 为 年 ,计算 时 取 m=1; 

A 是 “缩放 系数 ” ,表示 对 不 同 数据 源 进行 标准 化 处 理 , 使 得 时 间 信 息 箭 的 计算 结果 具有 可 比 性 (一 般 A 
值 可 以 取 为 基础 变化 单位 ;以 本 文采 用 的 NDVI 数据 为 例 , A 即 由 时 序 NDVI 数据 统计 得 到 的 基础 变化 单位 
A=0.02) 。 

通过 时 间 信 息 炉 可 以 反映 某 一 时 段 内 植被 覆盖 在 时 间 维 的 变化 强度 情况 : 矿 值 越 大 反映 出 这 一 段 时 间 内 

其 变化 强度 越 大 ,反之 则 表示 变化 强度 越 小 。 分 析 时 间 信 息 灼 计算 方法 ,由 于 计算 时 对 遥感 观测 数据 进行 
了 排序 处 理 ,4 吉 果 只 能 反映 出 植被 覆盖 的 变化 强度 信息 ,而 无 法 体现 出 其 在 时 间 维 的 变化 趋势 信息 齐 基 于 此 ， 
本 文 提出 了 时 间 序 列 信息 箭 (time-series information entropy ) 来 反映 植被 覆盖 在 某 忆 时 间 范 围 内 的 变化 趋势 信 
息 。 以 研究 中 所 用 的 NDVI 数据 为 例 , 具 体 的 计算 见 式 (8): 


了 1 : 
H’ = SB (vim 一 时 > log, Re xm i . 
4,0 >0 
式 中 ,«, 表示 像 元 第 年 的 NDVI 值 ; sgn 表示 符号 函数 , sgn(9) = 410, 六 = 0 其 他 参数 合 义 同 上 。 
,0 民 R0 


万 值 为 正 表示 在 某 一 时 段 内 植被 覆盖 的 变化 呈 增 加 趋 殷 8 反 之 则 表示 植被 覆盖 呈 减 少 趋势 。H' 绝对 值 
越 大 ,表明 增加 (或 减少 ) 的 趋势 越 明 显 。 

利用 上 述 时 间 信 息 粮 .时 间 序 列 信息 粹 可 以 挖掘 植 钻 覆盖 在 时 间 维 的 变化 特征 ,定量 、 客 观 地 反映 出 长 时 
序 植被 覆盖 的 变化 强度 和 变化 趋势 信息 。 同 时 ,在 计算 得到 中 述 粹 值 的 基础 上 ,对 其 直方 图 进行 分 割 ,进而 得 
到 植被 覆盖 的 变化 等 级 分 布 图 ,可 为 区 域 生态 环境 的 保育 和 治理 提供 更 为 方便 宏观 的 参考 。 
2.3 线性 回归 分 析 法 

目前 的 研究 中 多 采用 线性 回归 分 析 法 来 反映 区 域 长 时 间 序 列 植被 覆盖 的 变化 情况 下 ,本 研究 将 新 建 
立 的 时 间 信 息 烂 及 时 间 序 列 信息 业 的 结 时 污 线 性 回归 分 析 方 法 对 比 ,以 分 析 新 方法 的 优势 。 

线性 回归 分 析 法 以 线性 拟 合 得 到 的 斜率 值 ( slope) 表达 植被 覆盖 的 变化 强度 和 变化 趋势 。 具 体 而 言 , 斜 
率 的 绝对 值 越 大 ,表示 植被 的 变化 强度 越 大 ,反之 则 变化 强度 越 小 ;斜率 值 为 正 ,表示 植被 覆盖 呈 增 加 趋势 , 反 
之 则 呈 减 少 趋势 。 翰 计算 方法 风 式 (9): 


(9) 


式 中 ,ww 为 监测 年 数 ; NDVL 为 任意 像 元 第 i 年 的 NDVI 1 值 。 
2.4 ”时 间 信 息 粹 方法 的 验证 比较 
为 验证 本 研究 所 提 新 方法 的 科学 性 和 有 效 性 ,本 文选 取 延 河流 域 为 试验 区 ,将 此 方法 与 目前 长 时 序 植被 

盖 变 化 检测 中 最 常用 的 线性 回归 分 析 法 进行 比较 。 长 时 间 序 列 植被 覆盖 的 变化 情况 与 区 域 的 生态 环境 等 
密切 相关 。 随 着 人 类 活动 .气候 变化 等 综合 因素 的 影响 ,区 域 植被 覆盖 的 变化 情况 越 来 越 复杂 ;不仅 有 变化 趋 
势 比 较 明显 的 区 域 , 还 有 变化 趋势 虽然 不 显著 但 在 时 间 维 的 变化 波动 性 较 大 的 区 域 。 对 于 长 时 间 序 列 植被 覆 
盖 的 变化 检测 而 言 ,其 重点 内 容 是 客观 .科学 地 反映 出 研究 区 某 段 时 间 内 的 变化 强度 和 变化 趋势 特征 ,进而 为 
区 域 的 生态 环境 建设 和 管理 提供 辅助 决策 信息 。 基 于 此 ,本 文 将 以 长 时 间 序 列 植被 覆盖 的 变化 强度 和 变化 趋 
势 检 测 为 标准 ,选取 代表 性 样本 点 ( 表 1 、 图 6) ,来 综合 对 比 回归 分 析 法 与 本 文 提出 算法 的 计算 结果 ,验证 和 
说 明 时 间 信 息 粹 方法 的 客观 性 和 科学 性 (详细 的 内 容 见 3.3 节 ) 。 
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表 1 代表 性 样本 点 的 实际 数据 


Table 1 Data of representative sample points 


年 份 样本 点 1 观测 值 样本 点 2 观测 值 样本 点 3 观测 值 样本 点 4 观测 值 
Year Samplel Sample2 Sample3 Sample4 
2000 0.8269 0.3862 0.4734 0.4835 
2001 0.8335 0.4126 0.4276 0.5009 
2002 0.8369 0.5188 0.4154 0.5556 
2003 0.8423 0.4157 0.4111 0.5112 
2004 0.8436 0.5169 0.4246 0.5196 
2005 0.8412 0.5072 0.4382 0.4825 
2006 0.8387 0.5065 0.3373 0:5623 
2007 0.8395 0.6021 0.3415 045015 
2008 0.8426 0.5687 0.359 08969 
2009 0.8433 0.6798 0.2956 0.5551 
2010 0.8516 0.5907 0.3154 0.5643 


3 ”研究 结果 与 讨论 


3.1 线性 回归 分 析 法 的 结 

利用 线性 回归 分 析 法 计算 得 到 2000 一 2010 年 延 河流 域 植被 覆盖 的 变化 情况 。 结 果 发 现 ,此 时 间 范 围 内 
研究 区 的 植被 纱 盖 总 体 呈 上 升 趋势 (图 2) 。 流 域 北 部 和 东南 部 植被 履 凑 的 上 升 趋势 较 大 ,西南 部 的 上 升 趋势 
较 小 。 植 被 覆盖 呈 下 降 趋 势 的 区 域 主要 分 布 在 流域 的 中 部 、 西 北部 , 且 中 部 植被 的 下 降 趋势 较 大 。 

利用 线性 回归 分 析 法 的 判定 系数 RR 来 检验 所 得 斜率 值 的 可 信 程度 。 计 算 结果 发 现 , 判 定 系数 较 大 的 区 
域 主要 分 布 在 延 河流 域 的 北部 .东南 部 ,表明 这 些 区 域 NDVT 回归 分 析 斜 率 值 有 较 高 的 可 信 度 (图 3)。 判 定 
系数 较 小 的 区 域 主要 分 布 在 流域 的 中 部 西南 部 西北 部 ,表明 这 些 区 域 的 计算 结果 可 信和 度 较 低 ,并 不 能 客观 
地 反映 植被 覆盖 的 变化 情况 。 


图 例 

本 0.035 一 0.072 
四 | 0.020 一 0.035 
|] 0.010 一 0.020 
四 0 一 0.010 

国 -0.025 一 0 
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i <0.15 
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36°30" 
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109°00’ 109°30’ 110°00" 110°30'E 109°00’ 109°30’ 110°00" 110°30 禄 
2 延 河流 域 NDVI 回归 分 析 斜 率 值 分 布 3” 延 河流 域 NDVI 回归 分 析 判 定 系 数 分 布 图 


Fig.2 Distribution of slope values in Yanhe watershed Fig.3 Distribution of determination coefficient in Yanhe watershed 


3.2 ”时 间 信 息 信 方 法 的 结 

利用 时 间 信 息 箭 方法 ,计算 得 到 2000 一 2010 年 延 河流 域 植被 覆盖 的 时 间 信 息 粹 .时 间 序 列 信息 粹 的 分 布 
情况 ,分 别 见 图 4、 图 5。 

结果 表明 ,此 时 间 范 围 内 流域 的 植被 覆盖 总 体 变化 强度 较 大 (图 4)。 延 河流 域 的 土地 利用 类 型 多 为 耕 
地 草地 :2 , 近 十 年 来 由 于 “退耕 还 林 、 还 草 ” 等 政策 的 实施 ,该 区 域 的 植被 覆盖 整体 变化 较 大 。 而 流域 的 南 
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部 和 西南 部 为 林地 ,此 时 间 范 围 内 植被 变化 强度 较 弱 , 故 其 时 间 信 息 烂 值 也 较 小 。 本 文 研究 结果 与 延 河流 域 
的 实际 情况 相符 合 。 


37°00N 
37°00'N 


36°30" 
36°30" 


1 1 1 1 
109°00’ 109°30’ 110?00' 110°30 禄 109°00’ 109°30" 110°00" 110°30'E 


图 4 ” 延 河 流域 植被 覆盖 的 时 间 信 息 粹 分 布 图 5 ” 延 河 流域 植被 覆盖 的 时 间 序 列 信息 炳 分 布 图 
Fig.4 Distribution of temporal information entropy in Yanhe Fig.S Distribution ofstime-series information entropy in Yanhe 
watershed watershed 


研究 表明 ,2000 一 2010 年 间 植 被 覆盖 呈 增 加 趋势 的 区 域 主要 分 布 在 流 焉 的 北部 .东北 部 和 东南 部 ,反映 
出 这 些 区 域 的 植被 覆盖 情况 变 好 (图 5) 。 而 流域 的 中 部 、 西 背部 和 西南 部 的 植被 覆盖 主要 呈 减 少 趋势 ,需要 
重视 此 类 区 域 植被 覆盖 的 保育 工作 ,避免 进一步 恶化 。 
3.3 算法 对 比 


以 长 时 间 序 列 植被 覆盖 的 变化 强度 和 变化 趋势 检 了 a 
测 为 标准 ,来 综合 对 比 回归 分 析 法 与 本 文 提出 算 潜 的 计 -| 二 司 0010 二 0020 
算 结果 ,验证 和 说 明 时 间 信 息 炉 方法 的 客观 性 得 科学 S| Se 是 os 
性 。 研 究 依据 延 河流 域 植被 覆盖 的 斜率 值 和 判定 系数 “| ro 


> 


的 分 布 情况 ,分 别 从 变化 趋势 平缓 < 上 升 和 是 降 区 域 中 ， 
选取 4 个 代表 性 样本 点 的 NDVI 数据 (图 \6) 。 

对 于 上 述 样 本 点 的 实际 数据 《 表 1) ,分 别 利用 线性 
回归 分 析 法 和 本 文 提出 的 时 间 往 息 人 方 法 计算 ,所 得 结 


36°30" 


果 如 表 2 所 示 。 109°00’ 109°30’ 110°00’ 110°30 福 
对 于 前 3 个 释 本 点 的 实际 数据 ,可 以 容易 地 看 出 ， 0 
as 、 i 6 性 观测 点 位 置 示 意图 
2000_2010 第 同样 本 乱 1 的 NDVI 数据 变化 非常 平稳 ， enh 


Fig.6 Location of typical observation points 


样本 点 2 增加 的 趋势 明显 ,样本 点 3 减少 的 趋势 显著 。 
相应 的 ,从 表 2 中 我 们 可 以 得 到 :由 回归 分 析 法 得 到 的 斜率 值 其 判定 系数 较 大 。 故 利用 线性 回归 分 析 法 得 到 
的 斜率 值 有 较 高 的 可 信和 度 ,能 够 反映 出 这 段 时 间 内 的 植被 覆盖 变化 特征 。 同 样 的 ,本 文 提 出 的 时 间 序 列 信息 
炉 志 反映 出 这 些 变化 趋势 :时 间 序 列 信息 焙 的 绝对 值 越 大 ,表示 变化 趋势 越 明 显 , 且 时 间 序 列 信息 焙 为 正 , 表 
示 植 被 覆盖 呈 增 加 趋势 ;反之 ,表示 呈 减 少 趋势 。 

对 于 样本 点 1.4, 由 回归 分 析 法 得 到 的 斜率 值 相 同 , 但 其 变化 的 波动 性 却 有 着 很 大 的 差别 。 对 于 植被 覆 
盖 的 此 种 变化 情况 ,线性 回归 分 析 法 难以 客观 地 反映 出 时 间 维 的 变化 特征 。 以 样本 点 4 数据 为 例 , 利 用 该 方 
法 得 到 的 斜率 值 的 判定 系数 非常 小 (0.106) ,也 说 明 其 可 信 度 极 小 ,不 能 客观 地 表征 植被 的 变化 特征 。 随 着 人 
类 活动 的 加 剧 、 气 候 变 化 等 综合 因素 的 影响 ,导致 区 域 植被 覆盖 的 变化 波动 较 大 、 变 化 日 趋 复杂 。 以 延 河流 域 
为 例 ,统计 图 3 中 所 示 的 判定 系数 分 布 情况 ,可 知 其 中 46.9% 的 区 域 由 回归 分 析 法 得 到 的 斜率 值 其 判定 系数 
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较 小 (小 于 0.65) , 故 利用 回归 分 析 法 得 到 的 斜率 值 难以 客观 地 反映 出 此 类 区 域 的 植被 覆盖 的 变化 规律 。 利 
用 本 文 提 出 的 时 间 信 息 烂 方法 能 够 很 好 地 表征 出 此 类 变化 特点 :样本 点 4 的 时 间 信 息 和 值 大 ,反映 出 这 段 时 
间 内 NDVI 变化 的 波动 性 较 大 ;而 样本 点 1 的 时 间 信 息 和 值 具有 0.3557 ,表明 其 变化 的 波动 性 非常 小 。 综 上 ， 
我 们 能 够 得 到 :本 文 提出 的 时 间 信 息 箭 能 够 客观 表征 植被 覆盖 变化 的 波动 性 特征 。 


表 2 计算 结果 比较 


Table 2 Comparison of the results 


算法 样本 点 1 计算 值 样本 点 2 计算 值 样本 点 3 计算 值 样本 点 4 计算 值 
Method Results of samplel Results of sample2 Results of sample3 Results of satnple4 
斜率 值 Slope values 0.0015 0.0237 -0.0156 0.0015 
判定 系数 Determination coefficient 0.6558 0.7673 0.8015 0i0106 

时 间 信 息 炉 Temporal information entropy 0.3557 3.7075 3.0973 2.5343 

时 间 序 列 信息 炉 


0.8204 2.0994 =1,5576 1.3599 


Time-series information entropy 


综合 以 上 分 析 ,虽然 线性 回归 分 析 法 可 以 简单 .直观 地 反映 2000 一 2010 年 延 河 流域 植 被 覆盖 的 变化 情 
况 。 然 而 ,由 于 该 方法 的 前 提 假 设 是 植被 在 一 段 时 间 范 围 内 呈 线 性 变化 特 诛 ,对 于 植被 变化 波动 较 大 的 区 域 ， 
线性 回归 分 析 法 却 很 难 客观 .科学 地 表征 这 段 时 间 内 植被 覆盖 的 变化 特征 绚 耐 本文 提 出 的 时 间 信 息 丧 方法 表 
征 长 时 间 序 列 植被 覆盖 变化 信息 有 着 独特 的 优势 ,能 够 客观 地 反映 费 研 究 区 植被 覆盖 的 变化 强度 和 变化 趋势 
特征 。 

更 进一步 的 ,通过 对 粹 值 分 布 的 直方 图 进行 分 割 , 得 到 植被 窗 盖 变化 等 级 分 布 图 ,可 为 区 域 生态 环境 的 保 
护 和 治理 提供 更 为 方便 、 宏 观 的 参考 。 以 延 河流 域 为 例 和 中 反映 植被 覆盖 变化 强度 信息 的 时 间 信 息 粹 值 分 布 
图 进行 图 像 分 割 得 到 阔 值 A=1.68( 来 确定 基本 不 变 的 医 域 闯 对 反映 植被 覆盖 变化 趋势 信息 的 时 间 序 列 信 息 
烂 进行 分 割 得 到 阔 值 B=1.96( 判定 增加 的 程度 是 否 明 显 ) .C=-0.73( 判 定 减少 的 程度 是 否 严重 ) 来 进一步 反 
映 研 究 区 植被 覆盖 的 变化 等 级 。 综 合 以 上 的 内 容 $2000 一 2010 年 延 河 流域 植被 覆盖 变化 情况 可 分 为 以 下 五 
种 等 级 。 

时 间 序 列 的 和 遥感 检测 信息 反映 近 Toi 年 流域 生态 环 
境 偏好 (图 7)。 统 计 结 果 表 明哲 项 得 到 2000 一 2010 
年 延 河流 域 植被 覆盖 的 变化 等 级 以 增加 为 主 , 占 流域 面 
积 的 80.7%; 明显 增加 的 区 城 占 流 工 面积 13.9% ,主要 
分 布 在 流域 的 东北 部 和 东南 部 ;植被 覆盖 减少 的 区 域 占 
比 2.4% ,主要 分 布 在 流域 的 西部 和 西北 部 ;严重 减少 的 
区 域 占 比 1.19%5 主 要 分 布 在 流域 的 中 部 和 西南 部 ,是 需 
要 重点 的 坐 态 恢复 与 治理 区 域 ;基本 不 变 的 区 域 占 流域 
面积 的 工 9% 。 
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图 7 ” 延 河流 域 植被 覆盖 变化 等 级 分 布 图 
本 文 利用 时 间 信 乱 烂 方法 计算 了 2000 一 2010 年 延 Fig.7 Distribution of vegetation change levels in Yanhe 


河流 域 的 植被 覆盖 变化 信息 ,厘清 了 其 时 空 变化 特征 。 ”watershed 
研究 表明 , 延 河流 域 的 植被 覆盖 情况 总 体 得 到 改善 , 反 
映 出 该 区 域 的 “退耕 还 林 、 还 草 ” 政 策 得 到 了 有 效 实施 ,与 延 河流 域 的 实际 情况 相符 合 。 然 而 ,流域 的 中 部 、 西 
部 的 植被 覆盖 呈现 恶化 趋势 ,应 重视 此 类 区 域 植 被 覆盖 的 保育 工作 ,加 强生 态 环 境 的 修复 和 治理 力度 。 

本 文 提出 的 时 间 信 息 炉 方法 与 回归 分 析 法 相 比 ,能 够 更 加 客观 地 表征 长 时 间 序 列 植被 禾 盖 的 变化 强度 和 
变化 趋势 信息 ,可 为 区 域 生 态 环境 的 保护 和 管理 提供 更 为 科学 的 决策 参考 。 研 究 成 果 可 为 长 时 间 序 列 遥 感 影 
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像 的 植被 覆盖 变化 检测 提供 新 方法 .新 思路 。 研 究 旨 在 揭示 某 一 时 间 范 围 内 区 域 植被 覆盖 的 整体 变化 特征 ， 
尚未 涉及 其 中 的 “突变 ”问题 ,未 来 研究 需要 利用 时 间 信 息 焙 方 法 ,以 期 揭示 植被 覆盖 时 空 变化 的 影响 机 制 。 
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